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Exponential derivatives and generalized reproduction factors 


A very short introduction 


Abstract. The exponential derivative of a differentiable path c: J > G in a Lie group acts through right 
or left multiplication in G and offers so an intrinsic description of the infinitesimal change, in contrast to the 
classical derivative é = dc/dt that acts in tangent space TG~,). An elementary example is the exponential 
derivative of a smooth strictly positive functions f :J 4 R+ into the abelian multiplicative group G = R+. 
As a prototype, this is covered in detail. 

Basic properties, similarities and differences to the classical derivative are outlined and calculation rules 
are listed. In the Appendix, an application to Mathematical Epidemiology is presented. 
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1. Introduction 


The exponential derivative is nothing new. It is even used extensively, just in these times of the Corona 
pandemic, without being explicitely mentioned as such. Instead we encounter the term reproduction 
factor almost every day in the news, and it tells us how many people are infected on average from 
one contagious person, cf. [yDr]. One might wonder what this seemingly discrete concept can have 
to do with a derivative. There is no smooth variable visible in its definition, and so a connection 
to a derivative seems puzzling at first glance. But time enters naturally with the transmission time, 
which is the average time 7;, during that an epidemic with reproduction factor R spreads from one 
infected person to R susceptible ones. The fraction u(t) of newly infected individuals per time unit, 
the so-called incidence, then increases or decreases with time according to 


u(t + At) & RAVE .a(e) (1.1) 
in higher than first order with At > 0, thus introducing a natural time scale for the epidemic spread. 


The mapping R/” = exp [x 4 /u(t)] is in this case the exponential derivative of the incidence u(t). 
As an analytic one-parameter subgroup, it approaches u at time t within the multiplicative group R4, 
rather than in the additive group of tangent space (isomorphic to R, in this case) like the classical 
differential. 

Relation (LJ) calls for generalization and Lie groups (beside quantum mechanics, cf. page 4) provide 
a natural framework, thanks to the existence of the exponential map. That is what this study will be 
mainly about. In essence, our task boils down to making the natural connections to the exponential 
map explicit and relating them to what is well known. However, since this is an introductory study, 
we also build a viable bridge to calculus by treating the elementary case G = R; in more detail. 

The exponential derivative expresses the infinitesimal change of a C°-path c:J 4 G in any Lie 
group G by inner multiplication, i.e. using a factor in G, in contrast to the classical differential which 
acts by addition in tangent space. Consistent with that, equation (1.1) can be read as a ‘finite quotient’ 
instead of a finite difference equation. 

In an elementary abelian model, for example based on G = RY‘, the exponential derivative may 
describe the evolution of a pandemic, a population growth, the energy freed in a nuclear reaction, or 
any other dynamical process that follows some kind of chain reaction. Typically, the dynamic equations 
exhibit transient exponental increase or decrease with a time-dependent base (as exemplified in the 
Appendix). 

Applications for the general Lie group framework remain to be found. In particular, we think of 
dynamical systems with multiple interacting components. 

Last but not least, the whole theory is so natural, clear, and almost self-explanatory in many steps, 
especially in its elementary setting (cf. Section 4), that one might wonder why it is not part of any 
advanced Calculus course. 
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2. Basic ideas 


The terms differentiable, smooth, or C° will be used synonymously with differentiable in every order. 
For instance, a differentiable path is a C™-path and a differentiable manifold is a C’°-manifold as 
defined in standard textbooks of Differential Geometry, e.g. [Ddn], [Ist], [Hel]. 


In classical Analysis, the differential of a C™-path c: J > R” int € Jisa linear map Dc;(x) : R > R” 
and in particular a homomorphism from the additive group of R into that of R” such that 


< e(s)|, = <[Dela) + c(t) (2.1) 
which implies that Dc;(x) + c(t) approximates c in t in higher than the first order, that is 
c(t + At) = De(At) + c(t) + O(At) . (2.2) 
In analogy to (LJ) this is written 
c(t + At) © De, (At) + c(t). (2.3) 


In these equations, the differential Dc; acts directly in R” thanks to the canonical identification of R” 
with its tangent space. The situation is obviously different in a differentiable manifold where no such 
natural identification exists, quite apart from the fact that there is usually no addition or another 
binary operation in an arbitrary C-manifold, with which equations (2.2) and (2.3) would make 
sense. Again, things are different, if the manifold is a Lie group G. Then the differential Dc; is still 
operating linearly in the tangent space TG,;) , however, via the exponential map exp : g > G, the 
differential Dc; gives rise to an analytic one-parameter subgroup in G, i.e. an analytic homomorphism 
Ec,(xz) :R— G into the connected component of the unit element 1g € G, such that Ec;(x) acts 
from the left (or alternatively from the right) directly in G and Ec;(At)-c(t) approximates c(t+ At) 
for At > 0 in higher than first order, viz. precisely 


d d 
< o(s) |, = L[Bee(2) - (tI (2.4) 
in correspondence to (2.1). In that case we write in analogy to (2.3) 
c(t + At) ® Ec;(At) - c(t). (2.5) 


Based on this formal correspondence, which is made precise in the following section, Ec;(x) is named 
the (left) exponential derivative of c at t. 


Notice that relation (2.4) extends naturally to C°-paths c:J > M ina differentiable manifold if the 
connected component of 1g of the Lie group operates diffeomorphically on M (or at least on an open 
environment of c(J) ). 


From now on we just write simply gh instead of g-h for the binary operation in G. 


3. Geometric analysis 


Given a Lie group, i.e. a differentiable manifold G with an analytic binary operation which makes 
G a group and given a differentiable path c:J— G defined on an interval J C R, we first state the 
following 


THEOREM 1 


1. For every t € J, there exists an unique analytic homomorphism Ec; : (R,+) 9G, «4 Ec (x), 
i.e. an analytic one-parameter subgroup of G, such that 


5 OO: = ae [ Ecx(x) c(t) |, - (T1-1) 


2. The mapping 


is differentiable. (T1-2) 
In particular, for every fixed « € R, the mapping Ec: t+ Ec;(a) is a differentiable path J > G. 
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Proof: Define the function 


c:JIxJI—G (3.1) 
(s,t) + c(s)c h(t). 
Then c is obviously differentiable as a composition of differentiable maps. Also, c(t,t) = 1g (the 
unit element of G'), hence e(s,t)|, = £c(s) ew), € TGi, = g and the exponential map 
Ea(z): RG d (3.2) 


ue> exp (x=-c(s, t)], ) 


is well defined for every ¢ € J and is an analytic homomorphism of the additive group of R into G. 
According to (3.2), Ec;(a) satisfies 


Fe(s)e"(O|, = ZEa()| (33) 


dx 
The latter is however equivalent to equation (T1-1), since for every g € G the right action of g on 
G is an analytic diffeomorphism, hence in particular the action of g = c(t) for every fixed t. The 
uniqueness of Ec;(x) then follows from the uniqueness of the exponential map exp(a v) , which as an 
analytic homomorphism of R into G satisfies a exp(x V)|o =v, for every v €g (cf. (Hell, chap II, 
Corollary 1.5). 


0° 


With @ denoting the Lie group of all differentiable paths c : J + G and Hom(R,G) the set of 
analytic homomorphisms from R into G, the relationship between Ec; and the exponential map is 
summarized in this commutative diagram 


Hom (R,G) 
en exp 
eCxJ 7 g= TGig 


GSO. 
(c,t) H £c(s) c(t)" ], 
Definition 1 


The unique function Ec: t+ Ec:(z) € Hom(R,G) satisfying equation of Theorem 1 is called 
the (left) exponential derivative of c in t (we feel somewhat tempted to name it the (left) fractional 
of c in t, in obvious analogy to the classical differential ...). 

Since we deal only with the left exponential derivative in the following, the adjective left will be 
omitted from now on. The ’right case’ is treated analogously. 


As usually, let the right and left translation and the conjugation with g € G be defined by 
Rg:G—G, Llg:G—G and o9:G—>G, respectively . 
hr> hg he> gh hr o9(h) :=ghg* 


As is well known, R, and L, are analytic diffeomorphisms and $9(h), also simply written A, is an 
analytic endomorphism of G. With these objects, the following derivation rules hold 
THEOREM 2 (BASIC RULES) . 


Let g €G and c,c~:J —G be differentiable paths. 
The following is easily veryfied: 


1. If c = const, then Ec = 1g (ie. the trivial homomorphism R52 4 1g ) ( ) 
2. For every nE€N, Ec” = (Ec)” . ( ) 
3. If c, c~ commute, then Ecc™ = EcEc™ . (T2-3) 
4. ER,c = Ec (T2-4) 
5. ELge = EcI (= Ed%(c) ). ( ) 
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It is obvious that, if c : ] —0o ,oo [> G isa differentiable one-parameter group, i.e. c(r+s) = c(r) c(s) 
for every r,s © R, then Ec, = c for every (ER. 


Corollary to (T2-3): If G is abelian, then the mapping c(t) 4 Ec,(1) is an endomorphism of the 
group of differentiable paths c:J ~G. 


Proof : We prove statement (T2-3), the remainder is left to the reader. 
Let t € J. Then, since c, c~ commute, 


Ecc? = exp{ az [e(s)e~(s)e~(t)* e(t) "JI, } 


d Mirsky ASE 
= exp{x— [e(s) c(t)" e*(s)e~() J], }- (3.4) 
Now we apply the following lemma, from which everything else follows trivially. 
Lemma. Let c1,c2:J—>G be differentiable paths and t € J, such that c(t) = co(t) = 1q@. Then 


d d d d 
eels 1 = 2 als = sete + elt , (3.5) 


Note that (3.5) still applies to non-commuting paths cy, co. In fact, it is easily seen that the differ- 
entiable mappings 


f:GxGHoG, g:J7~GxG, 9: I7AGxG, go: JA GxG, 
(u,v) HR uv s+ (c1(s) , c2(s) ) s+ (c1(s),1¢) s+ (le, ca(s)) 
satisfy fog=aam, fom =a and fog =e. 
Therefore i F . j 
els a df | 5) Pree oa df | 46,16) =, = df | 46,16) (350 t? ole) 
d d d d 
= Of 35 [Ceres 0) ag (0, qelt)! = lo. qe + f | got) qe Plt 
d d 
atl + ale: (3.6) 
With ci(s) := c(s)c(t)~+ and co(s) := c~(s)c~(t)~+ equation (3.4) becomes 
7 d 
Eccy = op (ag meals) 
= ola ge el +25 oe) 
= exp(x—cy|,) exp(x—ca|,) [since c, and cy commute] 
ds —'t ds “lt 
= Ec. Ect . (3.7) 


Remarks 


Since the exponential map of the additive Lie group R” is just the identity of R” , the exponential 
derivative of a path in R”, characterised by Theorem 1 (T1-1), coincedes with the classical derivative. 
This, besides the analogous properties listed in Theorem 2, also justifies the naming ‘derivative’ in 
the generalized sense of Definition 1. 

Apart from Differential Geometry, one finds equations of the form (T1-1) and hence some kind 
of ”exponential derivative” also in Quantum Mechanics, where it describes the infinitesimal change 
of a time-dependent evolution operator U(t), t € [0, oo[ within the group of unitary transformations 
U(H) of a separable Hilbert space H. 

In fact, for any suitably smooth family tt» H, of self-adjoint energy operators with a dense 
common domain DCH (ie. Hy is a time-dependent Hamiltonian), every state vector ¢(q,t) = 
U(t) d(q,0) € D evolves in H according to Schrédinger’s equation 

_O 

which yields, in analogy to (T1-1), 
d d @ 
— [U(s) o(a, 0) ]], = — [exp (=Hz) U(t) oa, 0) ]I, 
ds j ae a 

: x 
that is = (SI, - <, lexp (= Hi) Ul] 

d 


or also 4 Hy = 4U(s)U(t)~! l, weakly on D, in formal correspondance to (3.3). 


» weaklyonD, (3.9) 
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4. The elementary case G = Ry 


Let T € R;U{+co} and J be the interval [0, T'[. Also, let C be the convex cone of strictly positive 
functions in C°(J). Trivially, € is an abelian Lie group under function multiplication and the constant 
functions in @ form an abelian subgroup Const} C C. Therefore the factor group C; := C/Const; 
is well-defined and is via the map f +> f/f(0) isomorphic to the convex subgroup of @ in which all 
functions have initial value equal to one. We can thus identify C, with { f € C| f(0) =1} and denote 
by 71: f 4 f/f(0) the pertinent projection. In the same manner, C™(J)/Consty is identified with 
the subgroup C§°(J) := {f € C™(J) | f(0) = 0} of the additive group C™(J), the pertinent projection 
being 70: f > f — f(0). 

Let f € C andt € J. By Theorem [J] the unique function Ef;(x) : R > Ry satisfying 
is the exponential derivative of f in t,ie. Efi(z) = exp[xf(t)'4 f(s)|,] = exp[z + log f(s), ] 
and by the Corollary to (T2-3), the mapping 

BE:€—-¢€ 
f -— Ef; = Ef,(1) 


is an endomorphism of €. In the following we continue to write simply Ef; for Ef,(1) and agree 
upon these 
Definitions 
1. The unique operator J : C — (,, f ->o Jf, such that (EI) 
EoJ:€ >C =Ide:frof (ie. the identity of C), 
is called the exponential integral on €. 


Notice this implies that E’ is surjective and J is injective. 


2. Let a,b €J, fee, and J: C > C, as in (EI). (Elab) 
The exponential integral of f over the interval [a,b] is defined as 


Jf \a= Ife/ Ifa (= If()/ FF(@) )- 
3. Finally, under the same conditions as in (Elab) and additionally a < b , we call (EM) 


1 


Flay = (IF la es 


the exponential mean of f over [a, b] . 


The next theorem affirms that the just defined objects exist and that they are unique, i.e. they are 
well-defined. 


THEOREM 3 ( EXISTENCE AND UNIQUENESS OF Jf ETC.) 
1. There exists a unique exponential integral J : € —> C; as defined in (EI) 


and that is 
t 
Jf, = exp | log f(r) dr (T3-1) 
0 
2. The exponential integral over [a, b] CJ equals 
b 
ifs exp | log f(r) dr . (T3-2) 
3. Also, the exponential mean of f over [a, b] is 
_ 1 b 
fia] = exp 7 _ | log f(r) dr . (T3-3) 
Proof: Trivially, the composition of 
E:e€—->e é with J:C—>C; t 
ft) okEf= exp G, log f(t) f(t) 0 Jfp= exp | log f(r)dr 
0 


assigns to every f € C 


d t 
EoJ(f)= exp |, log exp | log f(t)dr =f. 


Hence Eo J = Ide. Similarly it is seen that Jo E | C, = Ide, . (4.1) 
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Furthermore J is an exponential integral as defined in (EI) and J is unique. In fact, by (T2-2) 
E:foE£Ef is homomorhic, ie. E(fg) = (Ef) - (Eg), for every f,g € ©. Moreover, E [| Ci is 
injective due to (4.1), as is also J: C > @,. It follows that for every f € C;, Ef =1 if and only if 
f =1, and that for any J; and Jz as in definition (EI) and g € C 
Eodig = Eodog — > Eodig/EoJgg = 1 

— (E 1 C1)(Jig/Jog) = 1 

<= Nig/ Jog = 1, 

—= Sig = Jag, 


which can only apply to every g € © if Jy = Jo. 
The remainder is just as straightforward. 


Note that identities (£1) imply that E [ C,:C; > @ and J : @C > C are bijective mappings, i.e. 
they are isomorphisms. Obviously this is not the case for E : € + C, since Ef = E(af), for every 
a €R,, in particular Ef = E(f/f(0)). 


The formal correspondance between the classical and the exponential derivative is already seen in 
some simple examples: 


D(const) = 0 
Daf + Bg) = aDf + BDg 


in particular: 


D(a+ f)=Df 


E(const) = 1 
E(f* - 9°) = (Ef)* - (Eg)? 


E(a-f)=Ef (ae€R,) 


Further properties are listed in the next 


THEOREM 4 (DERIVATION RULES) 


Let f,g € © and as above D denote 4. 
These rules are easily saa! 


1 E(f+9) = (Ef)T - (Eg) 
2. Be) -= Gays? 
3. E(fog) = [(Ef)og]?9 if 9) CJ, 
4 BE) = [64] if Df Ee, 
f(b = 
5. — = (fay), for [a,b] c [0, T[. 


Proof: Mere retrospection of the definitions and Theorem[3] We show the first and fifth statements 
and leave the rest to the reader. 


Ad 1. : 
E(f+g) = expDlog(f+g) = exp “E78 
a exp "exp 7 = (exp) 74s « (exp 22) 745 
= (Eft. (Eg) Te 
Ad 5. : 
f(b) mf (b) _ eeu )  JEf(b) ee _ 
f(a) 7 mf (a) 7 Emf (a) - JE a _ JEf la - (Ef (a,j ) 
(§@) (}Emf=Ef) (fElab) ($EM) 


In the Appendix, these concepts are applied to earthly reality in a mathematical model of epidemic 
spread. 
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5. Appendix 


A ’real world’ problem 


In the model of epidemic spread presented here, we consider strictly positive C'°°-functions on a 
time interval J:= [0,T[, T € Ry U{+oo}. Remember these functions form a convex cone and a 
multiplicative group C C C™(J) and C; = { f € C| f(0) =1} is a convex subgroup of €. Also, it 
is easily seen that 


EoJ=Ide and JoE}@; = Ide, [ cf. (4.1) ], (A.1) 


where FE: € +€ and J: C€ > €1, respectively, are the mappings 


t 
Ef, = exp £108 f(s) i and Jf, = exp [ log f(s) ds. (A.2) 


Remember also that the exponential integral and the exponential mean of f over [a, b] C J equal 


b b 
Jf (®= exp | log f(s)ds and fia,b} = xD y ! / log f(s)ds, respectively. (A.3) 
a a 
For simplicity, the herd size is normalized to one in the following. The non-susceptible (’immune’) 
fraction of the herd is represented by a function h : J > [0,1] (0: no immune members, 1: total 
herd immunity ). More precisely, h(t) is defined as the unsusceptible fraction, comprising the immune 
members, but also those who are currently sick, and those who have died, i.e. h(t) = 1 — s(t), s(t) 
denoting the susceptible fraction. (Note that in our model, h(t) = 1 — s(t) is not equal to the sum 
i(t)+r(t) of the infected and recovered fractions, as is the case for example in the SIR model, cf. ; 
since we do not require that recovered members remain immune.) 

Finally, the incidence, defined as the fraction of symptomatic (’acute’) new infections per time, is 
represented by the function u : J > Ry. 

The reproduction factor indicates how many members are infected on average by one contagious 
member. By itself, it says thus nothing about the underlying time scale. Yet time enters naturally with 
the transmission time T;, which is the average time during that one sick member infects R susceptible 
ones. Then obviously u(t+7;,) = R-u(t) and we can state the following 


Proposition 1 (Reproduction factors) 


With Eu;z as above, the reproduction factor, usually defined as R := weed , equals: 


R — ( Bupet+Ti] Ver, (A.4) 
ie. RT is the exponential mean of Eu over the transmission interval [t, t + Ti, ]. 


Proof : Immediate consequence of Theorem 3 (T3-3). 
Remarks 


1. Since Bu : t ++ Fu, is continuous, there always exists a time T € [¢,¢+Ty,], such that 
Eu, = Eupt+r,,] , and hence, due to (A.4), 


Ri = Bu, (A.5) 


2. Trivially, the transmission time T;, introduces a natural time scale in which T;, is the unit of 
time. Hence in this scale, for a certain 7 € [t,t +1], 


R= Eu, . (A.6) 


From now on we work in this natural time scale, i.e. with T;, = 1. Moreover, we identify R = R; 
with Hu; = exp 4 log u(t), which amounts to refining the definition of R by replacing the exponential 
mean of Fu over a transmission cycle by its current value at time t. This exactly is summarized in 
the following 


Proposition 2 . In the natural time scale (in which T;, = 1), the reproduction factor is the expo- 
nential derivative of the incidence. 


8 STEFFEN HEIN 


Every mathematical model makes simplifications. An essential assumption of ours is a homogeneous 
environment. This means, to put it simply, that the pathogene finds the same conditions for spread- 
ing throughout the herd. In particular, singular spreading events or new infections from outside are 
neglected. As far as these conditions are met, the essence of the disease spread is captured by the 
following model equations, which relate the immune and fatal fractions to the incidence 


t t 
h(t) = h(0) + w; | exp ( )u(r) dr + wy | u(r) dr , (A.7) 
0 0 
C 

100 — A 100 — A™ 
The first integral represents the recovered uninfectable (’steril’) fraction of the herd and the second 
integral stands for the members that died up to time t. The decay time constant T; limits the mean 
duration of immunity after infection, and the weighting factors w; and wy introduce the percentage 
A of the asymptomatic (’silent’) cases, B of steril-immunized cases, and C of fatal cases into the 
equation (so if all infected survivors become non-infectiously immune, then obviously B + C = 100). 

Under the assumed homogeneous conditions, the exponential derivative of the incidence is directly 
proportional to the susceptibility. By virtue of proposition 2, this means 


R, = Eu, = s(t) Huo = Ro[1 — A(t)] 


= Ro [1 h(0) —w, | exp (>) u(r) dr 2 w, | Natale (A.8) 


T—t 


4 


with w; = andwy = 


and so the incidence becomes, in virtue of (AI), (EI), and (A.8) 
t 
u(t) = u(0) J Eu, = u(0) exp | log Eu, dr 
Sch (AL), 
“i 0 : a 
= u(0) exp | log { Ro [1 — h(0) — w; | exp ( )u(0) dd — wv; [ u(0)d0]} dr, (A.9) 
0 0 0 


where u(0), h(0) and Ro represent the initial incidence, the initial herd infection, and the reproduction 
factor at time t = 0, respectively. 

Equation (A.9) describes the development of an epidemic under homogeneous conditions. Other 
elements, such as a vaccination rate, are easily introduced into equations (A.7-9). In the AMADEUS 
program this and other features are implemented. 

The exorbitant dynamics that an epidemic can develop at higher reproduction factors is also 
reflected in the derivative of u(t), which follows directly from (A.9) 

t t 

hia] = u(t) log { Ro [1 — h(0) — wi f exp é = Saute dt — wv, [ u(r) dr]}. (A.10) 

dt '* 0 T; 0 
We have integrated equation (A.9) for basic reproduction factors varied from Ro = 0.5 to Rp = 3.0 
using a numerical routine implemented in the AMADEUS program, cf. . The initial conditions are 
h(0) = 0 and u(0) = 10~° (in common parlance: ‘incidence one’, i.e. one sick person in a population of 
100.000), and we assume T; = oo and B+ C = 100 percent, that is every surviving infected member 
ends up permanently immune. We also allow A = 20 percent asymptomatic cases. 


If 


4 
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Fig. 1. The incidence 


60 


The graphics give an idea of what is happening if an epidemic with higher reproduction factor is 
let loose, i.e. if it is not contained by suitable measures. For Ro = 2.7, and always under the above 
assumptions, a significant portion of the population would be simultaneously sick at the peak. 
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Fig. 2. The herd infection 


Under the above conditions, with Ro = 2.7 initially and no further action taken, everyone would 
eventually go through with the infection, be it alive or dead. The entire population would end up 
immune (the cuts in the two upper curves indicate when 100% herd immunity is attained). Assuming 
a mortality rate of 0.2%, for example, in a population of 82 million, as in Germany, more than 160,000 
people would die. In such conditions one can only hope for any effective vaccination, otherwise the 
reproduction factor would have to be kept low with painful containment measures. 


6. Conclusions and outlook 


The exponential derivative of a differentiable path c: J G in a Lie group has been studied and it 
was demonstrated how it acts within the group through right or left multiplication and in this way 
provides an intrinsic description of the infinitesimal change, in contrast to the classical differential 
that acts in tangent space. 

Basic properties have been outlined and calculation rules derived. The elementary case of the 
multiplicative group G = R, together with the smooth strictly positive functions f : J — R; has 
been treated in more detail, and applied to a numerical model of epidemic spread. 

Applications of the general Lie group framework should to be further explored. We are particularly 
thinking about chemical reactions or dynamical systems with multiple interacting components, but 
this must be left to future research. 
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